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We propose a method for systematically finding ground states of spinor Bose-Einstein condensates 
by utilizing symmetry properties of the system. By this method, we can find not only an inert state, 
whose symmetry is maximal in the manifold under consideration, but also a non-inert state, which 
has lower symmetry and depends on the parameters in the Hamiltonian. We establish the symmetry- 
classification method for the spin-1, 2 and 3 cases at zero magnetic field, and find a new phase in 
the last case. Properties of vortices in the spin-3 system are also discussed. 
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I. INTRODUCTION 

Classification of ordered states based on symmetries 
has been employed in many areas of physics, chem- 
istry, and mathematics. In ordered states, symmetries of 
the system at high temperatures are spontaneously bro- 
ken [l| . There exist several phases in quantum condensed 
systems with internal degrees of freedom, such as uncon- 
ventional superconductors, superfluid Helium three, and 
Bose-Einstein condensates (BECs) with spin degrees of 
freedom. The last ones are referred to as spinor BECs. 
Symmetries of the system are not completely broken 
in these phases, and once we know which symmetry is 
broken in the ground-state phase, we can immediately 
find what types of topological excitations, such as vor- 
tices, monopoles, and skyrmions, can be hosted in that 
phase [l-Ill- In this paper, we discuss how to find the 
ground state of a BEC with internal degrees of freedom, 
from a point of view of symmetry classification. 

Here, we briefly explain the concept of the symmetry- 
classification method for the case of spinor BECs. In the 
mean-field approximation, we assume that all atoms are 
Bose-Einstein condensed in a single-particle state, that is 
the order parameter of the system. For a spin-i^ system, 
the order parameter is a {2F + l)-component complex 
spinor: xjj = [ipp t4'f-1t ■ ' ' where T denotes 

the transpose. Then, the ground-state order parameter is 
obtained by minimizing the mean-field energy functional 
E[xIj] with respect to Vmi i-^-j 



Slpm 



= 0, m = F,F-l, 



-F. 



(1) 



The above set of equations gives the multi-component 
Gross-Pitaevskii equation for a stationary state. In gen- 
eral, to find the ground state, we must solve a set of 
2F -\- 1 nonlinear coupled equations. Though this pro- 
cedure works well when F is small [5|49 || , the calculation 
becomes very involved for large F |10l . . 

One reason for the complexity of the problem is that 
there are an infinite number of solutions to Eq. ^ asso- 
ciated with symmetry breaking. In general, a system 
under consideration has a certain symmetry, which is 



spontaneously broken in the ordered phase. For the case 
of spinor gases, the Hamiltonian is invariant under the 
global gauge transformation, global spin rotation, and 
time reversal. The mean-field energy is also invariant 
under these transformations on xjj. If we find a solution 
to Eq. ID), the order parameters obtained by applying the 
gauge transformation, spin rotation, and/or time reversal 
to the solution also satisfy Eq. ([ij . A set of order param- 
eters obtained by such transformations is called an orbit. 
For example, the direction of the spontaneous magnetiza- 
tion in a ferromagnetic state is arbitrary in the absence of 
an external field. All ferromagnetic states having differ- 
ent directions of magnetization belong to the same orbit 
and should be identified as the same class of states. 

A set of operations under which the order parame- 
ter remains invariant constitutes an isotropy group. The 
isotropy group characterizes the symmetry of the state, 
and its conjugacy class provides a convenient label to 
classify the individual states according to their symme- 
tries. Moreover, such a symmetry consideration gives fur- 
ther clues for finding the ground-state order parameters. 
According to Michel the gradient of the energy 

functional with respect to the order parameter vanishes 
in the direction along which the order parameter changes 
its symmetry. It follows that if we find a stationary solu- 
tion by restricting the order parameter space so that the 
order parameter has a certain symmetry, the obtained 
state is always stationary in the whole order parameter 
space. This theorem greatly simplifies the procedure for 
finding a stationary state in comparison with direct solu- 
tion of Eq. ([T])- In particular, in some cases, there is only 
one solution (orbit) which has a certain symmetry. Such 
a state, which is called an inert state, is always stationary 
and robust against a change in interaction parameters. 
Inert states have been obtained from the symmetry con- 
sideration in p- and d-wave superconductors Il3l . 11411 , su- 
perfluid helium three [l^, and spinor BECs [l6l Il7j. On 
the other hand, a non-inert state depends on the param- 
eters in the interaction energy, and energy minimization 
must be invoked to find it. 

In this paper, we discuss the symmetry-classification 
method and apply it to spinor BECs with spin-1, 2, and 
3 bosons at zero magnetic field. For the cases of spin-1 
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and 2 BECs, the obtained results agree with those found 
in the previous works For the case of a spin-3 

BEC, however, the systematic method enables us to find 
a new phase, which exists in a very narrow region in the 
space of the scattering lengths and has eluded the previ- 
ous works [lol . ITsI l . In the spin-3 system, there are many 
phases that have discrete symmetries, leading to various 
kinds of vortices as in the case of a half-quantum vor- 
tex in the spin-1 polar phase (l9j and a 1/3- vortex in 
the spin-2 cyclic phase [lO, [2lj. We discuss properties 
of vortices in spin-3 BECs and show that the quanti- 
zation unit of the mass circulation depends on the in- 
teraction parameters in some phases. The study on a 
spin-3 BEC has been motivated by the experimental re- 
alization of BEC of ^^Cr atoms [H, US]. Examples in- 
clude the phase diagrams in the presence or absence of 
an external field [ifl, [lli, those under a light- induced 
quadratic Zeeman energy [23 | or conserved magnetiza- 
tion j25j, and phase separation under an external mag- 
netic field [1^ . Possible vortices in each phase are inves- 
tigated in Ref. [l^ ■ Recently, spinor properties of spin-3 
^^Cr BEC have been observed |27l. [28|. 

This paper is organized as follows. In Sec. |TT1 we 
describe a method for finding a stationary point of an 
arbitrary function on a smooth manifold, and estab- 
lish mathematical notations used in the present pa- 
per. In Sec. mil we describe a general procedure of 
the symmetry-classification method in spinor BECs. In 
Sec. IIVI we carry out this procedure for spin-1, 2, and 
3 BECs. In particular, for the case of a spin-3 BEC, we 
point out a new phase which has eluded Refs. [l3, E^- In 
Sec. IIVI we discuss properties of vortices in spin-3 BECs. 
In Sec.|Vl we make concluding remarks. In the appendix, 
we explore stationary states with discrete symmetries in 
spin-3 BECs. 



II. SYMMETRY-CLASSIFICATION METHOD 

Our symmetry-classification method is based on the 
following Michel's theorems dlHlli]. 



A. Michel's Theorem 

We consider a real smooth function / on a smooth 
manifold A4. Let G be a group of operations which au- 
toniorphically map A4 to itself and do not change the 
value of /: 

G={ge AutM\f{gx) = f[x) for e X}, (2) 

where kxiiM. denotes the group of automorphisms on 
M.. For spinor BECs, / is the mean- field energy, M is 
the order-parameter manifold, and G is a group of gauge 
transformations, spin rotations, and time reversal (see 

sec irmi) . 



An orhit M.o{x) of a; S is defined as the trajectory 
of a point x on the manifold under G: 

Mo{x)^{gxeM\geG] ^M. (3) 

By assumption, / takes on the same value on all the 
points in JV[q[x). If a; is a stationary point in /, Aio{x) 
is also stationary. What we need to find is not a station- 
ary point X, but a stationary orbit M.o{x)- An isotropy 
group Hx is a set of operations that do not change x: 

= {g e G\gx = x} c G. (4) 

It is clear that Hx is a subgroup of G. It can also be 
shown that the isotropy groups of points on the same 
orbit are conjugate to each other: 

Hgx = {g' e G\g'gx = gx} 
= {g' e G\g~^g'gx ^ x} 

= {9g"g-' e G\g" e H^} 

= 9Hxg-'. (5) 

Here, H and H' , which are subgroups of G, are conjugate 
to each other if and only if there exists g G G such that 
H' = gHg~^. If two points on different orbits share 
the same isotropy group, the orbits of such two points 
are considered to be of the same type, and we classify 
the types of orbits according to the conjugacy classes of 
subgroups of G. In other words, for each conjugacy class 
of a subgroup of G, we obtain a set of orbits. Such a 
union of orbits is called a stratum A4s{x)', x and x' € Ai 
belong to the same stratum, if and only if their isotropy 
groups arc conjugate to each other. Clearly, Aio{x) C 
Ms{x) C M. 

Embedding the manifold Ai in an n-dimcnsional Eu- 
clidean space M" , where n > dimTM , the gradient of / at 
X = (xi, X2, • • • , Xn) is defined as 

y? ft ^ fdfdf df\ 

Here, VmUx) describes the direction of the steepest- 
ascent vector which is tangent to the manifold Ai . Michel 
has proved that '^Mf{x) is tangent to the stratum 
A4s{x), i.e., the gradient of / vanishes in the direction 
along which the symmetry of the state changes 0, HI] . 
Moreover, since / is a G-invariant function, VMf{x) is 
zero in the direction of Mo{x). Hence, we obtain the 
following theorems: 

Theorem 1 (inert state). If an orbit is isolated in the 
stratum, the orbit is stationary. 

Theorem 2 (non-inert state). If an orbit is not iso- 
lated in the stratum, we define a submanifold A4h C A4 
such that 

Mh = {x e M\hx = x br'^h e H}, (7) 

where _ff is a subgroup of G that characterizes the stra- 
tum under consideration. Let fn be a real function which 
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is the same as / but whose domain is restricted on J\4h- 
Then, the stationary point of on Mh is always a 
stationary point of / on A^. 

It follows from Theorem 1 that all G-invariant functions 
/ on a manifold A4 have a common stationary orbit. The 
corresponding state is called an inert state. Theorem 2 
is instrumental in finding non-inert states. 



B. Procedure 

Following the above two theorems, the procedure to 
find a minimum of a G-invariant function / on a manifold 
M is summarized as follows: 

1. Classify ah subgroups of G according to conjugacy 
classes. 

2. Let H be an element of a conjugacy class. Find such 
X G Ai that it is invariant under H. 

3. (inert state) If x is uniquely determined, then x is 
a stationary point of /, and the corresponding orbit 
Aioix) is a stationary orbit. 

4. (non-inert state) If x is not uniquely determined, then 
calculate the minimum of / in the submanifold Mh = 
{x S Ai\hx = X for ^/i G H}. A stationary point 
X € Mh also a stationary point of / in the whole 
space of M . 

5. Finally, compare the values of / for the obtained sta- 
tionary states and find the lowest one. 

We emphasize that the above procedure works well for 
the case in which stationary states have a certain sym- 
metry. If H = 1, i.e., if no symmetry remains, the above 
procedure amounts to solving Eq. ^ directly. We there- 
fore do not consider the case oi H = 1. In the absence of 
an external field, all ground states of spinor BECs with 
spin F ~ 1,2 and 3, superfluid '^He, and p- and rf-wave 
superconductors have remaining symmetries. 



III. GENERAL PROCEDURES FOR 
SYMMETRY CLASSIFICATION OF SPINOR 
CONDENSATES 

In this section, we describe general procedures for ap- 
plying the symmetry-classification method described in 
the preceding section to the case of spinor BECs. 



A. Mean- field energy 

The mean-field energy of a uniform system of spin-i^ 
atoms with mass M at zero magnetic field is given by 

m 

+ i JJdrdr' ^ VZy-'ir-r') 

rmir)rm'ir')^Pn'ir')Mr)- (8) 

We consider only a short-range interaction for V™^V^ (r — 
r'), and ignore the magnetic dipole-dipole interaction. 
Then, the interaction potential conserves the total spin 
of two colliding atoms and can be approximated with the 
delta function as 

V:T'{r-r') = S{r-r')CZV', (9) 

2F .2 

^""'= E -^as{Fm,Fm'\Vs\Fn,Fn'), (10) 

S— 0,cvcn 

where as is the s-wave scattering length of the total spin 
S channel, and Vs ~ J2ms \S, Ms){S, Ms\ projects a 
pair of atoms onto the total spin S state. From the fact 
that the inter-atomic interaction is elastic and conserves 
the total number and total spin of particles, the Hamilto- 
nian is invariant under the global U(l) gauge transforma- 
tion, the S'0(3) rotation in spin space, and time reversal 
e ;30]. Hence, 

G = C/(1)0 X S'0(3)f X e (11) 

is the full symmetry of spinor gases. Here, the subscripts 
4> and F denote the gauge and spin symmetry, respec- 
tively. If the scattering lengths satisfy special relations, 
G can be enlarged. For example, if all as's are equal, 
G = SU {2F -I- 1) X 8. In this paper, however, we do not 
consider such exceptions. 

In the absence of a trapping potential or a long- 
range interaction, the ground state is uniform with 
fixed density n. Introducing a normalized spinor ^ = 
(Cf,Cf-i,-- - ,C-fV as Vm = V"Cm, thc grouud state 
is obtained by minimizing 

'^'[C] = "^^2^ ^ E ^n"™ CmCm'Cn'Cnj (12) 
mnni'n' 

subject to thc normalization condition (^^(^ = 1, where Q 
is the volume of the system. Using relations Y^g = 1 
and Es['5(5' + 1) - 2F{F + l)]Vs = 2F • F [J, Eq. (USD 
for = 1, 2 and 3 can be rewritten as 

S^''>[C]=co + c,\f\\ (13) 
£(^)[C] =co + ci|/|2 + c2|^oo|', (14) 

=Co + Ci|/|2 + C2|Aoo|'+C3El^2A/|', (15) 

M 
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respectively, where 

/ = ^CJ^mnC«, (16) 

mn 

Am = ^(0, 0|Fm, Fn)UCn, (17) 

'rnn 

Aim - ^(2, M|Fm, ^^n)C™C«, (18) 

mn 

are the magnetization per particle, the spin-singlet pair 
amplitude, and the spin-quintet pair amplitude, respec- 
tively, with F = {Fx,Fy,Fz) being the vector of spin-F 
matrices and {S, M\Fm, Fn) the Clebsch-Gordan coeffi- 
cient. The coupling constants co,ci,C2 and C3 are given 
in terms of the scattering lengths as (6l-[T]| 



(19) 



F = 


1 


Co 


= 2a2- 


1- flo. 








Cl 


= 02 - 


ao. 




F = 


2 


Co 


= 3a4 - 


h4a2. 








Cl 


= Qi — 


a2, 








C2 


= 804 - 


- 1002 4 


7ao, 


F = 


3 


Co 


= 2a6 - 


h 9a4, 








Cl 


= 06- 


04, 








C2 


= lOae 


- 21a4 


+ llao 






C3 


= 7a6 - 


- 1804 -+ 


lla2. 



(20) 



(21) 



where c, = (4F - l)Mci/(47r?i2). 

For F = 3, using the identity J2si^(^ + 1) " 2F(F + 
l)]2pg = (2F • F)2, the last term on the right-hand side 
of Eq. ([T5|) can be rewritten as 



2M 



A/ 



1 1 

7 " 18' 



1 
126 



^ ^ ^ ^ Cm {F^LFi')mnC,n 



(22) 



By decomposing F^^F^ in the last term into the symmet- 
ric and antisymmetric parts: 



F^F, 



^e^.A^A, (23) 



Eq. ^2$ reduces to 



1 



1 



M 



(24) 



where A/" is the spin ncmatic tensor defined by [K 



,r _Y^/-* {Ff^F^ + F^Ff, 



(25) 



By definition, is a real symmetric tensor whose trace is 
given by Tr AA = F(F -1-1). Then, the mean-field energy 
for F = 3 [Eq. ([T5|) ] is rewritten as 

f^'^IC] =5o + £i|/|' + S2|AooP + 53TrAA2, (26) 

where Ci's are related to Ci's as 

1 - 5 



Co = Co - yC3, Cl = Cl - — C3, 
52-C2-^C3, 53 = Y^C3. 



(27) 



B. Procedure for the case of spinor BECs 



Following the procedure in Sec. Ill Bl we first classify all 
subgroups of G given by Eq. ([TT|) according to conjugacy 
classes. The representation of each component of G in 
the spin-F manifold is given by 



[/(l)^ = {e**l|</)eR}, 

e = {i,r}. 



SO{i)F = {e-*^^"e-*^''''e-'^^^|a,/3,7 e M}, (29) 



(28) 
(29) 
(30) 



where 1 is the (2F -I- 1) x (2F + 1) identity matrix, a, /3 
and 7 are Euler angles, and the time-reversal operator T 
acts on ^ as 



(rc)m = (-i)"c 



(31) 



Since the eigenvalue of an arbitrary element g of 5*0(3) 
is given in the form of e^^ with A G M, the corresponding 
eigenstate is invariant under e~'^^g S t/(l)0 x S'0(3)_f, 
that is, it is invariant under a spin-gauge coupled op- 
eration. If there are two eigenstates that have different 
eigenvalues, these two states differ in the spin-gauge sym- 
metry. Hence, the procedure 1 and 2 in Sec. IIIBI are 
rephrased as follows: 

1'. List all subgroups of SO(?>)f- 

2'. Let Fl' be a subgroup of 50(3)i?. Find simultaneous 
eigenstates C of all elements of each H' . 

For procedure 1', it is known that SO{2) is the only 
continuous subgroup of 50(3), and that the discrete sub- 
groups of 50(3) are given as follows [3l| : 

Gn- The cyclic group of rotations about a symmetry 
axis through angle 2'Kk/n with fc = l,2,-- - ,n — 1. 
The group is isomorphic to Z„. 

Dn'. The dihedral group generated by the elements of 
Gn and an additional rotation through tt about an 
orthogonal axis. 

T: The point group of the tetrahedron composed of 4 
three-fold axes and 3 two-fold axes. 

O: The point group of the octahedron composed of 
3 four-fold axes, 4 three-fold axes, and 6 two-fold 
axes. 
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subgroup 


generators 


SO{2) 




Cn 




Dn 


Cnz 5 C2x 


T 


Csz, C^y^x + z 


O 


Ciz, C2,x + z 


Y 


Csz, C'2,2a;+(l+\/5)z 



TABLE I. Generators of the discrete subgroups of SO{3), 
where Cn,ui^x+uiyy+ui^z denotes a rotation about the direc- 
tion u: — (uJx,iOy,uJz) through 2ii/n as given by Eq. (|32|) . For 
example, generator 2a;+(i+\/5)z ™ icosahedral group de- 
notes a TT-rotation about u) = (2, 0, 1 -f %/5). 



Y: The point group of the icosahedron composed of 6 
five-fold axes, 10 three-fold axes, and 15 two-fold 
axes. 

Since all rotations through a fixed angle about different 
axes are conjugate to each other, we choose a represen- 
tative element H' in a conjugacy class of each subgroup 
of 5*0(3)^ so that the highest symmetry axis is parallel 
to the z axis. The generators of such representative el- 
ements of the subgroups of SO{Z)f are summarized in 



denotes a rotation about 



Table m where Cn, 

the direction u> = {uJx,LOy,i-Uz) through 27r/?i: 



Cn^Lj^x-l-ujyy+LJzZ — exp 



.ujxFx + oJyFy + uj^F;, 27r 



(32) 



For example, the matrix elements of Cm and C2x hi a 
spin-i<" system are given by 



[(7 1 , — ~i2TTm/ns: 
\C2x\ra7n' — ( 1) ^ra. — ni'- 



(33) 
(34) 



Next, we calculate simultaneous eigenstates of all gen- 
erators of each subgroup H' . First, we consider the 
case of H' = 5*0(2). The eigenstate of the generator 
Fz is uniquely determined as [C^™°^]m = 5mrao- Here, 
we neglect the overall phase factor, since e**^^^™") G 
■A^o(C^'"°^) for ^0 G M. We also disregard negative mo, 
since the order parameter is obtained by applying 

the time-reversal or the spin-rotation operator to ^("'o); 



'(-mo) 



e'^^C2xC 



(mo) 



(mo) 



For mo 7^ 0, the isotropy group C G of C*""^ 

given by 



H 



(mo) 



{gimoog- 



i[2mo7+(-F+mo)7r] jiT-y-j ^gg-j 



where 



^2,(cos^)x+(siii^)y — ^ ^2x^ 



(36) 



is a 77 rotation about an axis in the x-y plane, and a and 
7 are arbitrary real numbers. Here, we have used the 



following relations: 



gimoag-iFjQ^(mo) 



^(mo)^ (37) 



i-rrF 2imo7/'(™o) 



gi[2mo7+(-F+mo)7r] |--y-^(mo)-|^^ ^gg-j 



Since and its time-reversal 7~C are transformed 

to each other through a spin rotation, the time-reversal 
symmetry is broken in these states. On the other hand, 
the mo = state has the time-reversal symmetry, which 
is decoupled from spin rotations in the isotropy group: 



H 



(0) 



(39) 



Note that the mp = state also has the Z2 symme- 
try [loj : the order parameter is invariant under the tt 
rotation about an axis perpendicular to the z axis, and 
the isotropy group can be written as TJ'"' = D^o x 0, 
where Doc denotes the dihedral group of order n = 00. 

For the case of discrete subgroups, we shall calculate 
simultaneous eigenstates and the corresponding isotropy 
groups for each subgroup H' in Sec. |lVl Most of them 
are not determined uniquely, and we minimize the energy 
in the submanifold spanned by the simultaneous 
eigenstate for each H' . We note that the obtained state 
for a given H' might have higher symmetry than H' . For 
example, since C„ C Dn, the stationary point in the sub- 
manifold with symmetry C„ may have the Dn symmetry. 
To identify the symmetry of the obtained state, it is con- 
venient to use the Majorana representation |32| . which 
we explain in the next subsection. 



C. Majorana representation 

Majorana invented a geometrical representation of a 
general spin-F state [11]. A state of the spin-F system 
can be specified by providing a symmetric configuration 
of 2F spin-1/2 systems, expect for an overall phase. Since 
the state of a spin-1/2 system can be described by the 
unit-sphere Bloch vector, the state of a spin-F system 
can be described by 2F vertices on the unit sphere. This 
representation helps us identify the symmetry structure 
of a spin-F condensate (33| . 

We consider a polynomial of degree 2F for a given 
order parameter ^: 



PPiz) 




(40) 



a=0 



Then, the 2F complex roots of P^^' (z) ~ give 2F ver- 
tices on the unit sphere through the stereographic map- 
ping z = tan(6'/2)e"^. For the cases of spin F = 1, 2 and 
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3; the polynomials P^^\z) arc given, respectively, by 



2. Continuous symmetry 



P^'\z)=C*z' + V2CoZ + Ci, (41) 

+ VTbCiZ^ + V6C2Z + Cs- (43) 

Using the Majorana representation, we can immedi- 
ately find some inert states For example, a spin-2 
BEC is described with 4 vertices. When these 4 ver- 
tices form a tetrahedron, the corresponding state has the 
tetrahedral symmetry. On the other hand, since an oc- 
tahedron and an icosahedron have 6 and 12 vertices, re- 
spectively, no spin-2 state has octahedral or icosahedral 
symmetry. The octahedral symmetry appears in systems 
with F > 3, and the icosahedron symmetry in systems 
with F >6. 

For the time-reversed state TC, the polynomial is given 
by 



2F 



(44) 



Q=0 



If z = zo is a root of (z) = 0, then z ~ —^/zq is a 

root of P^^ {z) ~ 0, which corresponds to the antipole of 
zq on the unit sphere. Hence, the time-reversed state is 
described with antipoles of vertices of the original state. 



IV. SPECIFIC EXAMPLES IN SPINOR 
BOSE-EINSTEIN CONDENSATES 

In this section, we apply the procedure discussed in the 
previous sections to spin _F = 1, 2, and 3 systems. 



A. Spin-1 

1. Mean-field energy 

A spin-1 spinor BEC is described with a three- 
component spinor C ~ (Cii Co, C-i)"^- The scaled mean- 
field energy for a given order parameter ^ is given by 



£W[C]=co + ci|/| 



(46) 



where 



f+ = f, + ify = V2CiCo + V2CoC-i, (47) 
/. = |CiP-|C-ir (48) 



There are two inert states that have continuous 
isotropy groups: 



F:C^'^ = il,0,Of, 
P:C'"'- (0,1,0)^, 



(49) 
(50) 



where the former is the ferromagnetic state, while the 
latter is the polar (or antiferromagnetic) state. The 
isotropy group of these states are given by substituting 
(F,mo) = (1, 1) and (1,0) to Eqs. ([35]) and ([39]), respec- 
tively, as 



P : i/(")={e 



, e'^'u]} X e. 



(51) 
(52) 



where a and 7 are arbitrary real numbers. Substituting 
(^(1-0) jj-^ Yjq. (|46p . the mean-field energies are obtained as 



F: fW[C(^)]=co + ci, 
P : £(i)[C 



(0)1 



Co- 



(53) 
(54) 



3. Discrete symmetry 



We first consider the eigenstate of Cm whose matrix 
representation in the spin-1 manifold is given by 



Cnz = Diag[e 



, l,e 



\2ixlr, 



(55) 



where Diag means the diagonal matrix. A nontrivial 
eigenstate of Cnz that has more than two non-zero com- 
ponents exists only for n = 2 and is given by 

C2 : C^^^) -(yT^,0,V^)^ (56) 

with eigenvalue —1, where < 77 < 1/2. Note that we 
can always choose Ci and C-i to be real and positive 
without loss of generality, since the phase factors of these 
components can be removed by a spin rotation about 
the z axis and a gauge transformation. The geometric 
structures of C*'^'', C^^'*, and ^('-^2) g^j-g ghown in Fig. [TJ 

For the case of H' = D2, the order parameter should 
be an eigenstate of C2x- 



C2X 




(57) 



resulting in 77 = 1/2. This state is nothing but the polar 
state since C''^^^'? = 1/2) = (l/^/2,0, 1/V2)^ is related 
to C^°^ by rotation: e^'^/^g-^-f^^'^/^C^") = C'^'K?/ = 1/2)- 
In other words, C,'^^^ and C^'^^H'? = 1/2) are on the same 
orbit. This fact can also be understood by comparing the 
geometric structures of these two states (Fig. [1]). For the 
case of H' = C2, substituting the order parameter ([55)1 
in Eq. (j46| . we obtain 



f(i)[C(^^)] = co-fci(l-2,7)^ 



(58) 
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The stationary point of S^^^^ is at 77 = 1/2, and the corre- 
sponding state has the same symmetry as the polar state. 
Hence, there are only two phases in a spin-1 BEC: the fer- 
romagnetic phase and the polar phase given by Eqs. 
and (|50|) . respectively. 




FIG. 1. (Color online) Majorana representation of (a) ferro- 
magnetic state C^^\ (b) polar state C'^'i ^^'^ (c) state 
[Eq. (|56p ] with C2 symmetry. The color of the points shows 
the value of cos 6 = (1 — \z\^)/{l + \z\'^) according to the color 
gauge. For the case of the ferromagnetic state, the polyno- 
mial P^^\z) = has two degenerate roots at the north pole, 
which is indicated with x2 in (a). In (c), two vertices are 
on the y-z plane on the same latitude. They are both at the 
north pole for 77 = and on the equator for r; = 1/2. In (c) 
the C^^^\v — f/2) state has the same symmetry as the polar 
state. 



4- Phase Diagram 



Comparing the energies of the polar and ferromagnetic 
phases, wc obtain the phase diagram of the spin-1 sys- 
tem as shown in Fig. [21 The physics of the phase dia- 
gram is quite simple: from the mean- field energy (|46p. 
I/I should vanish for ci > 0, whereas it becomes maxi- 
mal (i.e., I/I = 1) for C]< 0; the former is polar and the 
latter is ferromagnetic p, Q . 



1st order 



ferromagnetic 

•X2 




polar 








FIG. 2. (Color online) Phase diagram of a spin-1 BEC. The 
phase transition at ci = is first-order because the symmetry 
changes discontinuously at that point. 



B. Spin-2 

1. Mean-field energy 

A spin-2 spinor BEC is described with a five- 
component spinor ^ = (C2, Cij Co, C-ii The scaled 
mean-field energy for a given order parameter ^ is written 
as 



f(2)[C] =Co + Ci|/|2+C2|A 



00 1 



where 

/+ 
/. 

^00 



2C2*Ci + V6CK0 + VeCoC-i + 2C-1C-5 



2IC2 
1 



IC1P-IC-1P-2IC-2I 



-^(2C2C-2-2CiC-i+Co)- 



(59) 

(60) 
(61) 

(62) 



2. Continuous symmetry 

There are three inert states that have continuous 
isotropy groups: 



F2 : C*^' = (1,0,0,0,0)^, 
Fl : C*'' = (0,1,0,0,0)^, 
UN : C*"' = (0,0,1,0,0)^, 



(63) 
(64) 
(65) 



where F and U N stand for ferromagnetic and uniaxial- 
nematic js^ . [ssj . respectively. The isotropy groups of 
these states are obtained by substituting {F,mo) — 
(2, 2), (2, 1) and (2, 0) in Eqs. ([25]) and respectively: 



F2 : 


H(2) 


Fl : 




UN : 


Hio) 



(2) ^ {e^^<^e-'^''',e^'m]T}, 



^-iF,a TTt 



u]} X e. 



(66) 
(67) 
(68) 



where a and 7 are arbitrary real numbers. Substituting 
^(2,1,0) ggj)^ 

we obtain 



F2: £(2)[C(^)]=co+4ci, 
Fl: £:(^)[C(i)]=co+ci, 



UN: £:(2)[c(")]=co + ^ 
5 



(69) 
(70) 
(71) 



The Majorana representations of F2, Fl, and UN states 
are shown in Figs. [3Ia), [31[b) and[31Jc), respectively. 



Discrete symmetry 



Ci. 

Since a spin-2 system is described with four vertices 
in the Majorana representation, the symmetries that a 
spin-2 BEC may have are T, I?4,3,2, and C4_3^2- For each 
symmetry group, we seek stationary points of f'-^^IC]. 
The Majorana representation of the obtained states are 
shown in Fig. |3] (sec also Ref. [ssj). 
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(a) F2 

(1,0,0,0,0) 
'9X4 



(b) F1 

(0,1,0,0,0) 



(c) UN 

(0,0,1,0,0) 





(d)T (e)D4 

(1,0,0,i2,0)/i3 (1 ,0,0,0,1 )/J2 



(f)C4 

(JT^,0,0,0,Jri) 




(g)C3 

(Jl^,0,0,iry,0) 



G) D2 (k) D2 

(0<Ti<1/4) ;'°='(1/4<Ti<1) 



1^ 

cose 

1-1 



FIG. 3. (Color online) Majorana representation of spin-2 
BECs. (a) F2, (b) Fl, and (c) uniaxial nematic (UN) states. 
The F2 and Fl states have the SO{2) symmetry, whereas 
the symmetry group of the UN state is Doo ■ (d) Cyclic state 
with the tetrahedron symmetry. Four vertices in (d) form a 
regular tetrahedron, (e) D4 state with four vertices forming 
a square, (f) State with the C4 symmetry given in Eq. (|79p . 
Depending on j), the square moves between the north pole 
(rj = 0) and the equator (ry = 1/2). (g) State with the C3 
symmetry given in Eq. (|84|l . The triangle in (g) moves be- 
tween the north pole (77 = 0) and the south pole (r^ = 1), with 
one vertex fixed at the north pole, (h) and (i) represent the 
states (0, 1, 0, 1, 0)^/^2 and (1, 0, i\/2, 0, 1)^/2, respectively, 
(h) has the same symmetry as (e) , and (i) has the same sym- 
metry as (d). (j) C<°^> = W{l-il)/2, 0, ^, 0, ^{l-T^)/2f 
for Q < r) < 1/4 with four vertices forming a rectangle on the 
equator, (k) This state has the same order parameter as (j) 
but for 1/4 < 77 < 1. The state (^(^2) jjg^g ^j^g game symmetry 
as the D4 state when 77 = and 1/2, and as the UN state 
when 77 = 1/4 and 1. 



T: The generators of the tetrahedron group are C; 



3z 



spin-2 manifold are given by 

C3. = Diag[e^2V3^ g.4V3^ g.2V3^ g.4V3]^ 



c. 



exp — z- 



( 1 2^/2 2^/6 4\/2 4 \ 

2^f2 5 -2 -4^2 

2^/6 2\/3 -3 -2^3 2^/6 

4^/2 -2 -2\/3 5 -2^f2 

V 4 -4V2 2V6 -2\/2 1 / 



(72) 
(73) 

(74) 



The simultaneous eigenstate of these two operators 
is uniquely determined to be 



T: C'^^ = (1/^3,0,0,^273,0)^, 



(75) 



or its time reversal TCP^^ ~ 
(0, — ■\/3/2, 0, 0, l/-\/3)^, up to an overall gauge. 
The eigenvalues of ^^-^^ are e'^'^/'^ for C^z and 1 
for C2 ^x+z- This state is called the cyclic state. 
The isotropy group is generated by a set of the 
following operators: 

fjiT) ^ {e-^^-/^C,z,C,^^^^^,e-^^'m^T}. (76) 

where 7 is an arbitrary real number, and here and 
henceforth, we denote a set of generators of H by 
H. Though no spontaneous magnetization arises 
in the cyclic state, the time reversal symmetry is 
broken, because TC'"^'' 7^ C'"'^^ ■ One can also con- 
firm this fact from the Majorana representation as 
shown in Fig. [3jd): the antipoles of the vertices in 
Fig. El^d) form a time-reversed tetrahedron, but it 
does not coincide with the original one. The mean- 
field energy of the cyclic state is given by 



5(2) [^(T)] 



Co. 



(77) 



D4: A nontrivial eigenstate of 

C4. =Diag[-l,-z,l,i,-l] (78) 

is written without loss of generality as 

(v/r^,0,0,0,V^)^, (79) 

where < 77 < 1/2. When 77 — 1/2, this state 
becomes an eigenstate of C2X with eigenvalue 1. 
Hence, there is an inert state that has the -D4 sym- 
metry: 

Di : C^'""' = (1/v^, 0, 0, 0, l/V2f. (80) 

This state is often refereed to as the biaxial nematic 
state [13, UHl • The generators of the isotropy group 
and the mean-field energy for this state are given 

by 



f(^HC(^^)] = CO + I 



and C2 ^x+z whose matrix representation on the 



(81) 
(82) 
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respectively. Though the energy of is the 

same as that of the geometric structures of 
these states are different from each other as shown 
in Figs. He) and He). 

C4: The energy of the eigenstate of C4Z [Eq. ([7^ ] is 
calculated function of 77 as 

f(2) = CO +4ci(l- 277)2 + 1^,^(1 -,7). (83) 

5 

Equation (j83p has a stationary point at 77 = 1/2, 
resulting in the 1)4 state. The rj dependence of the 
geometric structure of the state given in Eq. (|79p is 
shown in Fig. [Sff). 

D3: The eigenstate of C^z is given by 

(71^,0,0,7^,0)^, (84) 

where < r; < 1 . The order parameter in the form 
of (0, y?7, 0, 0, \/l — rj)'^ is also a nontrivial eigen- 
state of C^z- However, this state belongs to the 
same orbit of that in Eq. ((84|) . because they arc 
transformed into each other by a spin rotation and 
a gauge transformation. There is no simultaneous 
eigenstate of C^z and C2x, i-e., there is no state 
that has the symmetry. 

C3: The energy of the eigenstate of C3 [Eq. (|84| ] is cal- 
culated as a function of 77 as 

=co + ci(2 - 377)2. (85) 

The stationary point of this function is 77 = 2/3, 
and the corresponding order parameter is identi- 
cal to the cyclic state [Eq. ([75)) ]. The geometric 
structure of the state given in Eq. ([M|) is shown in 
Fig.Hg). 

D2 ■ The matrix representation of C2Z is given by 

C2. =Diag[l,-l, 1,-1,1] (86) 
There are two simultaneous eigenstates of C'2z and 

C2x- 

Case (i): The order parameter 

(0,1/V2,0,1/V2,0)^ (87) 

is the simultaneous eigenstates of C22 and C2X with 
eigenvalues —1 and 1, respectively. However, this 
state has the same symmetry as that of the biaxial 
nematic (D4) state as shown in Fig. ^h). 

Case (ii): The other simultaneous eigenstate can 
be written as 

(/^,0,e-V^,0,/iI^)^ (88) 



where the eigenvalues of C2Z and C2X are both equal 
to 1. The mean- field energy of this state is given 
as a function of 77 and S as 

f^^) = Co + — [1 - 277 + 2772 + 2?7(1 - 77) cos 26] . (89) 
5 

Taking the partial derivatives of Eq. ([5^ with re- 
spect to S and 77, we obtain two stationary points at 
6 ~ 7r/2 and 77 ~ 1/2, and at S — and arbitrary rj. 
For the former case, the corresponding order pa- 
rameter, (1/2, 0, i/-\/2, 0, 1/2), has the tctrahedral 
symmetry as shown in Fig. ^i) . For the latter case, 
the order parameter is given by 

D2 : C^^^) = 0, VV, 0, ^i^) , (90) 

which has symmetry different from that of the other 
obtained state [Fig.^i)], as shown in Figs.Hj) and 
Hk). The energy for this state is calculated to be 

£(2)[C(^=)]=co + |, (91) 
5 

which does not depend on 77, implying that all 
the states described by Eq. (pi)) are degenerate. 
The uniaxial and biaxial nematic states are also in- 
cluded in Eq. ()90p , and they can be smoothly trans- 
formed to each other by changing 77 in Eq. ()90p . For 
a fixed rj, the generators of the isotropy group of 
^(^2) jg given by 

H^""'^ = {C2z,C2x,T}. (92) 

However, if we take into account the degrees of free- 
dom described by 77, the isotropy group of the state 
given in Eq. (HHj) is shown to be [Z2 x S'0(4)] x 
& [3^ , where x implies that the nontrivial element 
of Z2 does not commute with some elements of 
5*0(4). It has also been pointed out that the degen- 
eracy with respect to 77 is lifted if we take into ac- 
count quantum or thermal fluctuations [s^ . [35l [37| . 

C2' There are two nontrivial eigenstates of C22. 
Case (i): The order parameter 

(0,yT^,0,V^,0)^, (93) 

is the eigenstate of C22 with eigenvalue —1. The 
mean-field energy of this state is given by 

£(2) = Co + ci(l- 277)2 + 1^77(1 -77). (94) 

5 

This function has a stationary point at 7/ = 1/2. 
The corresponding state has the same symmetry 
as the biaxial nematic (-D4) state. 

Case (ii): The order parameter 

(a+,0,6,0,a-)^ (95) 
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is the eigenstate of C2z with eigenvalue 1, where a± 
and b are complex numbers. Here, we choose a± to 
be real numbers and rewrite these parameters as 

b = e^'^, (97) 

where 77,'^>0, ?7 + ^<l, and — tt < S < tt. The 
mean-field energy of this state is given by 

f (2) + 16ci^(l - ?7 - 

+ ^|l-2e-(l-e2'*)77p. (98) 

Taking the partial derivatives of Eq. (|M)) with re- 
spect to (5, 77, and ^, we obtain the same stationary 
solutions as those in case (ii) of the D2 symmetry. 

4- Phase diagram 

Comparing the energies of the obtained stationary so- 
lutions [Eqs. dM]), dlQl), dill), and wc obtain the 
phase diagram of a spin-2 BEC as shown in Fig. 01 In 
the region of nematic phase of Fig.|4l all states described 
by Eq. (|90p are degenerate, including uniaxial and biaxial 
nematic states. Our results agree well with those in the 
previous works Distinct from the case for = 1, 
the phase diagram for F = 2 is determined by the last 
two terms in Eq. ([5^ : Ci|/p and C2IA00P. Clearly, |/| 
can vary within < |/| < 2 for an F = 2 system. Note 
that 1^00 1 is proportional to the inner product of the or- 
der parameter and its time reversal: \Aqq\ = {TCyC/V^- 
It takes the maximum value of when the order pa- 
rameter has the time reversal symmetry, while it should 
vanish for the ferromagnetic state. Then, the ferromag- 
netic phase arises for ci < and C2 > 20ci, while the 
nematic state, which has the time reversal symmetry, be- 
comes the ground state for 02 < and ci > C2/2O. In 
the region of ci > and C2 > 0, the cyclic phase appears 
since both |/| and |Aoo| vanish in this phase. 

C. Spin-3 

1. Mean-field energy 

A spin-3 BEC is described with a seven-component 
order parameter C = (Cs, C2, Ci, Co, C-i, C-2, C-a)"^- The 
scaled mean- field energy is given by Eq. (|15p . Following 
the notations in Ref. [10[, we rewrite Eq. p5|) in terms 
of Boo = V7A00 and B2M = a/7^2M as 

2 

^C^\f\^ + Ca\Boo\^ + Cp l^2A/|', (99) 

M=-2 



ferromagnetic {F2) 

W 


cyclic (T) 


C2=20Cl / H 


nematic (D2) 





Ci 



FIG. 4. (Color online) Phase diagram of a spin-2 BEC. In 
the nematic (antiferromagnetic) phase, all states described 
by Eq. (|90p are degenerate. The phase boundaries on Ci = 
0, C2 =0 and C2 = ci /20 are all first-order because the sym- 
metry changes discontinuously across the boundaries. 

where — Ci, Ca = co/7, = C3/7, which correspond 
to 7, a, and /3 in Ref. [l3], respectively, 

/+ -V6C3*C2 + VT0C2C1 + 2\/3crco 

+ 2V3C0C-1 + VToc:iC-2 + Vec^aC-s, (loo) 

/. =3|C3p + 2IC2P + ICiP - IC-iP - 2IC-2P - 3|C-3P, 

(101) 



are the spin densities, and 



Boo 


=2C3C-3 - 2C2C-2 + 2CiC- 


^2 

-1 ^ Co: 


(102) 


B20 


=-^C3C^3 - V3ClC-l + 1 




(103) 


B2.±l 


= —^C±3Ct2 - ^C±2Ct1 


+ /fciiCo, 


(104) 



B2,±2 =y yCisCTi - \l yC±2Co + V2Cii. (105) 

2. Continuous symmetry 

There are four inert states that have continuous 
isotropy groups: 

C<^^ = (1,0,0,0,0,0,0), (106) 

C'') = (0,1,0,0,0,0,0), (107) 

C<^^ = (0,0, 1,0,0,0,0), (108) 

C(°) = (0,0,0,1,0,0,0). (109) 
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The isotropy groups of these states are given by substi- 
tutmg (F,mo) = (3, 3), (3, 2) and (3,1) in Eq. dSS]), and 
(F,mo) = (3,0) in Eq. ^ as 



as 



^(0) ^ {e-'-P^=",[/T'e"} X e, 



(110) 
(111) 
(112) 
(113) 



where a and 7 are arbitrary real numbers. Substituting 
^(3,2,1,0) (|99| . the mean-field energies are obtained 



f^(')[CW]=c, + 2c^, 



f(3)[^(0)] 



(114) 
(115) 
(116) 

(117) 



3. Discrete symmetry 



Since the state of the spin-3 BEC is represented by 
six vertices in the Majorana representation, the possible 
symmetries for a spin-3 BEC are O, T, -06.5,4,3,2, and 
(^6,5,4,3,2- For each symmetry group, we seek stationary 
points of f '^^^[C]. 

O: The generators of the octahedron symmetry are C^z 
and C2,x+z whose representations on a spin-3 man- 
ifold are given by 



Diag[i,-l,-i, l,i,-l. 



G 



2,x+z 



exp 



/ 


-1 


-x/6 




-2V5 


-^/T5 


-V6 


-1 \ 


-^/6 


-4 


-VTo 





yio 


4 


^/6 






-VTo 


1 




1 




-\/T5 














-2V3 





2V5 








1 


-2V3 


1 


VTo 


-VT5 




-^/6 


4 


-Vio 





Vio 


-4 




V 


-1 


V6 




2V5 


-VIE 


V6 


-1 / 



(118) 



respectively. The simultaneous eigenstatc of these 
operators is uniquely obtained as 



O : C^°^ = (0,l/%/2, 0,0,0, -l/\/2,0)^. (120) 



Here, the corresponding eigenvalues are equal to —1 
for both Caz and C2,x+z- Then, the generators of 



the isotropy group are given by 

H(0) = {e-C4., e"C2,.+., e^T}. (121) 
The mean-field energy is calculated as 

£(3)[^(0)j ^^^^ (122) 

T: The generators of the tetrahedron symmetry are 
Csz and C2 ^x+z whose matrix representations are 
given by 



12 



C: 



32 



Diag[l,e 



j:2ir/3 gi47r/3 ^1 



a 



exp 



1 

27 



/ 


-1 


-2^/3 


-2^15 


-4\/l0 


-4^15 


-8^3 


-8 \ 


-2V3 


-9 


-6^/5 


-2\/30 





12 


8^/3 




-2VT5 


-6V5 


-9 


2V^ 


12 





-4VT5 




-4VTo 




2\/6 


11 


-2\/6 


-2^/30 


4\/T0 




-4VT5 





12 


-2^/6 


-9 


6%/5 


-2^/15 




-8^3 


12 





-2\/30 


6^5 


-9 


2V3 


V 


-8 


8\/3 


-4VT5 


4VTo 


-2\/T5 


2\/3 


-1 / 



(123) 



These two operators have a unique simultaneous 
eigenstate 



A pa) 



B(D2; 



C(C2) 



(%/2/3, 0, 0, - V5/3, 0, 0, -\/2/3). 



(125) 



However, this state has the same symmetry as the 
O state. In fact, the Majorana representation of 
Eq. (|125p is obtained by rotating that the state in 
Eq. (fT20| 

For the symmetry groups of Z3„ and C„, we proceed 
in a manner similar to the case of the spin-2 BEC. We 
characterize the eigenstate of Cnz with a few parame- 
ters (77,^, etc.). For the Dq symmetry, the order param- 
eter is uniquely determined; therefore, this is an inert 
state. For the other cases, we rewrite £^'^^ in terms of 
the new parameters and find stationary points. The de- 
tailed calculations are described in Appendix El The re- 
sults are summarized in Table Ull in which we list all the 
obtained stationary states, together with their isotropy 
groups. We have obtained the analytical solutions for 
all stationary states, except for state C. For state C we 
have numerically calculated the energy by restricting the 
order parameter in the form of (a, 0, 6, 0, c, 0, d)^ with 
a, b,c,d G C. In the obtained state, a, b and d are real 
positive numbers and c is a real negative number. The 
geometric structure of the obtained states are shown in 
Fig. E] (see also Ref. [11]). 

Among the obtained stationary states, only states A, 
D, and Q possess the time-reversal symmetry. The time- 
reversal operator T is decoupled from spin rotations, 
which can also be understood from Fig. [SJ where the 
antipodal map does not change the configurations of ver- 
tices for the A, D, and Q states. The time-reversal op- 
eration changes the configurations of vertices for other 
states. In particular, the time-reversal symmetry is bro- 
ken in the B and E states, even though these states have 
no spontaneous magnetization, as in the case of the cyclic 
phase in a spin-2 BEC. Spontaneous magnetization arises 
in the states except for A, B, D, E, and Q (see Table. Hill) . 



4- Phase Diagram 

By comparing the energies, we obtain the phase dia- 
gram of spin-3 spinor BECs as shown in Fig. [51 Here, the 




□(O) 



E(D3) 



F (SO(2)) 




Q (SO(2)) 



R(C3) 



FIG. 5. (Color online) Majorana representation of stationary 
states of spin-3 BECs. The order parameter for each state is 
listed in Table HH 



energy of state C is calculated numerically. The phase 
diagram is almost consistent with that in Ref. [loj . How- 
ever, we have found a new phase J with the C4 symmetry, 
which has eluded the previous works [13, [ll] ■ We have 
also investigated the phase diagram by directly solving 
Eq. ([1]). and confirmed that no additional phase arises 
in the phase diagram. Because the method presented in 
the present paper can deal with only the states that pos- 
sess remaining symmetries, all the phases shown in Fig. [SI 
have certain remaining symmetries. 

Figure [HKa) shows the phase diagram for > 0. The 
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H' 


phase 


isotropy group H, H 


50(2) 


FF 


Eq. dllOl) 


50(2) 


F 
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50(2) 


P 


Eq. dm 




Q 


Eq. (ITT3ll 


O 


D 


{6^-04., e''^ (72,:.+., e"^' 


De 


A 


{e'''C;6.,e'"C2.,e"'r} 


Os 


H 





order parameter ^ 



04 



Di E 



O3 I 



O3 R 



D2 B 



O2 G 



C2 C 



{e''^''204.,02.r} 



(1,0,0,0,0,0,0) ^ 

(0,1,0,0,0,0,0) 

(0,0,1,0,0,0,0) 

(0,0,0,1,0,0,0) 

(0,1,0,0,0,-1,0)7^2 

(l,0,0,0,0,0,l)/^/2 

(y^,0,0,0,0,y^,0) 

eg 

^7 - 2(cg-3c^) 

(y^, 0,0,0, y^,0,0) 

_ 2cp 
^ 12c-y— eg 

(y^,0,0,yi±^,0,0, y^) 

48 



9c^ 

Cq + 3 



c +^ 



cg(72c.y— 25cg) 
12(3c-y —eg) 



e^(252c.^-25e^) 
12(12e.., —eg) 



Slca eg 



.ea+eg 



^g-i2,r/3(;.^^^ g-i./3y V6^| (0, 0, 0, 0, 0) 



{03.,e'"?72"/V} 
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12(3e^-eg) 



(l2C-y Cq, +16C-y — 9Cq ) 

12ci^ — c^c^ — CqC^ 



{02.,e"'^02.,e'"C4zr} 
{02.,e''^02.r} 



^7 2(3e^-eg) 

(a+,0,0,V^,0,0,a_) 
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9(12e^-e^eg-Caeg) 
^ 24e-^ Cq — 4c.y eg — 9ea eg 
^ ^ 18(12e:^-e^eg-eceg) 

m . / i-'? n i fl±E n . /I^ ^^^(isg-^-^seg) 

V^' V 4 I'-'i'y 2 ' ^' V 4 ' 3(3ea-4cg) 
eg 

7? = ^- 

' 4eg— 3cq 



V 



(0,a+,0,^,0,a_,0) 

(5cg — 6e.y ) (cq — 4e-y ) 
^ 8(eo eg — 4egc^+3e^) 
f _ eg(2ea-3e^) 
^ 4(e„eg — 4egc^ + 3c^) 

(a,0,6,0,c,0,d) 

a, 6, c, d G R, a, 6, d > 0, c < 



2Se„ e| -4e^c.^ (3e<j +25e^ ) + 12(3cq +4e^ )e^ 
12(eQCg-4e^cg+3ci^) 



numerically calculated 



TABLE II. Stationary states and mean-field energies for spin-3 BECs obtained by tlie symmetry-classification method. The 
second column indicates the phases discussed in Refs. [13, [3 except for phases J, P, Q, and R. State I in Ref. [l3| is refereed to 
as HH in Ref. P, Q, and R do not appear in the phase diagram of Fig. [6] The third column shows the isotropy groups. In 
A, D, and Q, the time-reversal operator T is decoupled from spin rotations in H, which implies that only these states have the 
time-reversal symmetry. In the calculation for Ds and D2, we have obtained the following states: (\/2, 0, 0, ±iy/E, 0, 0, \/2)"^/3, 
(0, ^,0,±yT0, 0,^,0)"^ /4, (^,0,^,0,^,0,^)^/4, (73,0,-^5,0,-75,0,^)^/4, and (1, 0, 7l5, 0, 715, 0, 1)^/4V^, 
which have symmetries of O, De, Dao, O, and De, respectively. For the D2 symmetry, we have also obtained the solution in the 
form of (a, 0, be'*, 0, fee"', 0, a)"^ {a,b,8 G R). However, this order parameter coincides with that of B by a 7r/2 rotation about 
the X axis. 



phase boundary between E and D is given by y = 5a;/3, 
while that between phase E and J is y = (252a; — 
25a;^)/[12(19x — 36)], where y = Ca/\c^\ and x = cp/\c^\. 
The B-C phase boundary is numerically obtained and 
well described by y = 3|a;|/(2|a;| -|-3). (The phase bound- 



ary y = 3|a:|/(3|a;| -I- 2) given in Ref. [T3l should read 
?/ = 3|x|/(2|a;|-|-3). Figure 1 in Ref. [l^ agrees with the 
latter one.) 

In the phase diagram for < [Fig. IHJb)], the phase 
boundary between H and G is y = (24x-(-5x^)/(364-24a;-|- 
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(a) Ci > 



6_ 

O 




12/5 
Cp/|Cy| 



(b)ci<0 





-3/2 



_t5_ 
O 



H, 1 

(degenerate) 


H 




G 


A 




\ 








D 



-6 --18/5 



Cp/|c,| 



FIG. 6. (Color online) Phase diagrams of spin-3 BECs for (a) 
Cj > and (b) < 0. The liglit-colored lines in (b) indicate 
second-order phase boundaries, while others show first-order 
phase boundaries. In (b), state H is degenerate with state I 
which can exist only when C/j/jc-^l < —6. 



5a;^), and that between A and FF is y = —25x/12 — 9. 
These results also agree with Ref. [10|. In the top left 
region of Fig. ^h) , states H and I, which have different 
symmetries, are degenerate. Distinct from the nematic 
phase in a spin-2 BEC [Eq. (|90|) ]. there is no intermedi- 
ate state between the H and I states, as pointed out in 
Ref. [Hi. 

The phase boundaries between G and A and be- 
tween FF and H are second-order because they can be 
transformed into each other by continuous changes in 
symmetry. Figure [7] shows the continuous symmetry 
change from A to G, and FF to H. In phase A, the 
order parameter has the Dg symmetry. On the A- 
G phase boundary, the order parameter of G becomes 
(0, ^/3,0,^/T0,0, v^,0)'^/4, where 77 = 5/8, and ^ = 0. 
This state has the Dq symmetry as shown in Fig. [7l[a 
As Co,/ 



vertices of FF lie at the north pole. As c^/|c^| decreases, 
one of the vertices remains at the north pole, and the 
other five move downwards while keeping the C5 symme- 
try [Fig.EJb)]. 



(a) A(D6) 

{0,i3,0,il0,0,i3,0)/4 

\ 



G (C2) 



M 




'if 




(b) FF(S0{2)) 
(1,0.0,0,0,0,0) 

X6 



H (C5) 



V 




FIG. 7. (Color online) Symmetry change at (a) A-G and (b) 
FF-H phase boundaries. In phase G, the four vertices move 
upward as indicated by the arrows and the symmetry in 
phase A breaks down. However, the six vertices are on the 
same plane and the C2 symmetry remains in phase G. In 
phase H, five vertices, which are on the north pole in phase 
FF, move downward while keeping the C5 symmetry. 

To discuss the underlying physics of the phase digram, 
we use the mean- field energy in the form of Eq. (|26p . 
[We used Eq. for the calculation, because the de- 
scription in terms of the order parameter is simpler for 
Eq. (HSl) than that for Eq. ^^.] In Fig. [3 we show the 
phase diagram in the parameter space of (ca/jci |, C2/|ci |) 
for (a) ci > and (b) ci < 0. For the case of F = 3, 
the phase diagram is determined by three terms: ci|/p, 
C2|^ooPj and c^TrN^. whose values for the obtained sta- 
tionary states are summarized in Table IIIII As we dis- 
cussed in Sec. IIV 1^41 |Aoo| takes its maximum when the 
order parameter has the time-reversal symmetry, that is, 
1^00 1 = l/v^ for A, D and Q. These states may be- 
come a ground state for C2 — > —00. They differ in the 
values of Tr N"^ . Since A/" is a 3 x 3 real symmetric ma- 
trix with trace F{F + 1) = 12, it has three real eigen- 
values, \i {i = 1,2,3), which satisfy J2i=i2 3^i ^ 1^. 



increases, the four vertices move upward as Then, Tr A/"^ = I]i=i,2,3 takes its minimum when all 



indicated by arrows in Fig. [Tl^a), and the Dq symmetry 
breaks down. However, the six vertices of G are still 
on the same plane and the C2 symmetry remains. On 
the other hand, on the FF-H phase boundary, the or- 
der parameter of H becomes (1,0,0,0,0,0,0)^ which is 
identical to FF. In the Majorana representation, the six 



three eigenvalues are the same, i.e., TrAf^ 



48 for 



(Ai,A2,A3) = (4,4,4). This is the case of the D phase. 
On the other hand, by noting that < Ai < = 9, the 
maximum value TrA/'^|max = 171/2 is achieved in the A 
phase in which the eigenvalues are given by (Ai, A2, A3) = 
(3/2, 3/2, 9). Hence, phases A and D arise for £3 < and 
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phase I /I 



TrA/'2 



(a)ci>0 H, I 

^ ' (degenerate) 



A 
B 
C 
D 
E 
FF 
F 
G 
H 



1 

|r/| 
numerically calculated 

1 

\r)\ 

3 

2 



4v/e(l-»]-C) 11-2^ 

T] 

I(HH) 77 

J 77 

P 1 

Q 1 



171/2 
6(9 + 277 



9(43- 677 + 277?2)/8 

171/2 

48 

24[2 - 477" + 577(1 - 0] 



R 



3(36 + 477 + 772)/2 
3(36 - 477 + r)^)/2 
[m-l2rj + l2if)/2 
123/2 

72 

6\/^(l -V-i) |1 - 277 - 2g| 9(19 - 30?7 + 2777')/2 



TABLE III. Magnetization j/j, singlet-pair amplitude |Aoo|, 
and nematic tensor Tr A/"^ for the obtained stationary states. 



£3 > 0, respectively, whereas phase Q does not appear in 
the phase diagram because its energy is always between 
those of A and D. Note that N^^^ — f^fu is the spin fluctu- 
ation, and hence, TiAf^ reflects the anisotropy of the spin 
fluctuation. The spin fluctuation is isotropic in phase A, 
while it is most anisotropic in phase D. 

In the region of C2/|ci| — ^ -l-oo and ca/jcil — >■ ~oo, the 
phase that has the minimum \Aqo\ and the maximum 
TtM^ becomes the ground state, i.e., the FF phase. In- 
terestingly, TrA/"^ becomes minimal for the F phase, al- 
though the magnetization is the second largest in this 
phase. Therefore, for the case of ci < 0, the F phase 
arises in the region of C2/|ci| — )■ -foo and ca/jcij — )■ -l-oo. 
On the other hand, there is no state that minimizes 
all I/I, |Aoo|, and TiAf^, simultaneously. Hence, many 
phases arise in the top right region of Fig. \^a.). In the 
limit of C2/IC1I +00 and C3/|ci| — > -foo, the J phase, 
which has I/I ^ 1/2, |Aoo| = 0, and TrW^ 48, be- 
comes the ground state. 



5. Vortices 

Vortices in spin-3 BECs have been classified in 
Refs. [l3,[l3- As explained in Refs. [T7l.[ll|. topologicallv 
stable vortices are classified in terms of the elements of 
the "lifted" isotropy group, which is a subgroup of the 
universal covering space of G, i.e., SU{2) x R x 0. Us- 
ing the symmetry-classification method, we can identify 
the isotropy group of the obtained state, and therefore, 
we can find what types of topological excitations can be 
hosted there. 

The order parameter far from a vortex core is generally 
described using gauge-transformation and spin-rotation 
operators as 



C(s) = 



rep: 



(126) 



where < s < 1 is a parameter describing a closed con- 
tour around a vortex, and ^j-cp is a characteristic order 




FF 


F 




D 


A 



0.0 \ 

2/15 



C3/IC1I 

FIG. 8. (Color online) Phase diagram in the parameter space 
of (c3/|ci|, C2/|ci|) for (a) £1 > and (b) ci < 0, obtained by 
transforming the parameters in Fig. [S] according to Eq. (|27p . 
In (a), the J phase exists in the region of large C2/|ci| and 
large C3/|ci|, as shown in the inset of (a). In the region on 
the right of FF in (a), H is degenerate with I which can exist 
only for C3/|ci| < —2/9. The curved phase boundaries in 
(a) are given by B-C: y = 168i(4 - Vox)/{2lx - 2), G-H; 
y = 252i(5i-2)/(45i^ + 12i-f4), andE-J: y = 567i/(4i-2), 
where x = C3/|ci| and y = C2/|ci|. The light-colored lines in 
(a) indicate second-order phase boundaries, while others show 
first-order phase boundaries. Phases A, D, and FF occupy the 
bottom left, bottom right, and top left regions, respectively, 
for both ci > and Ci < 0, because |/|, |Aoo|, and TrA/"^ 
take their minima or maxima in these phases (see Table ITlT)l . 



parameter for a state under consideration (see, for ex- 
ample, the order parameters in the fourth column of 
Tabic HH. For simplicity, we choose (j){0) = a{Q) = 
/3(0) = 7(0) = 0. Then, from the single- valuedness condi- 
tion for the order parameter [C(0) = C(l)]: the operator 
g#(i)g-^i=^,a(i)g-iF„/3(i)g-iF,7(i) must be an element of 

the isotropy group H. 

It is worth investigating the mass circulation of a con- 
densate, which is always quantized in a scalar BEG in 
units of h/M. In spinor BECs, mass circulation is not 
always quantized due to the contribution from the Berry 
phase caused by spin textures. The mass current of a 
spinor BEC is defined as 



E[C(VC™)-(VC)C™], (127) 



2Mi 



and its circulation is defined as a line integral of v„ 
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along a closed contour C: 

j> ■"mass ■ M. (128) 

In phases A, B, D, E, and Q, BECs have no magneti- 
zation, and mass current is proportional to the gradient 
of the overall phase 0: 

ft , ^ 

t^mass = -^V(/). (129) 

Moreover, these phases have the spin-gauge coupled Z2 
symmetry, namely, the order parameter is invariant un- 
der e'^'^C2x or e"C2^x+z- Then, the single- valuedness con- 
dition for these states requires that (/)(1) takes on an in- 
teger multiple of tt. It follows that k is quantized in units 
of h/{2M), which is one half of the conventional value. 
This situation is similar to the case of a half-quantum 
vortex in the spin-1 polar BEC [l^. Note, however, sev- 
eral topologically different vortices may have the same 
circulation. In particular, the A, B, D, and E phases 
host non-Abelian vortices because their isotropy groups 
are non-Abelian. Therefore, a topological charge that 
classifies each vortex in these phases is an operator (or a 
matrix) rather than a scalar quantity, as in the case of a 
1/3 vortex in the spin-2 cyclic BEC [g^,!!!!!!]. 

On the other hand, the other phases in Table Hill have 
spontaneous magnetizations, and the mass circulation is 
not simply quantized. For these phases, the mass current 
is calculated as 

t^mass = ]g - |/|(cos/3Va + V7)] , (130) 

where I/I is the amplitude of the spontaneous magne- 
tization given in the second column of Table. IIIIl and 
(sin /3 cos a, sin /? sin a, cos /3) describes the direction of /. 
Integrating Eq. (jl30p along a closed contour C, we obtain 
the following relation: 

- S{f) = 0(1) - |/|[a(l) + 7(1)] ^ 2nU, (131) 

where 

S{f) = \f\ jjl- cos p)Va-M (132) 

is the Berry phase due to a texture of /, and it is de- 
fined modulo 47r|/|. Since a and (3 specify the direction 
of the magnetization, they have to satisfy a(l) = 27rnQ 
and /3(1) = 0, where Ua is an integer. For the case of 
FF, F, and P, there is no discrete symmetry and the 
single-valuedness condition dictates that (/)(1) = 27™^ 
and 7(1) = 2Tm-y, where and are integers. Since 
I /I is an integer in these phases, /v defined in Eq. (|13ip 
should also be an integer. Due to the arbitrariness 
of the Berry phase, vortices are classified by integers 
mod 2|/|. For G, which has spin-gauge decoupled dis- 
crete symmetry, the single-valuedness condition leads to 



(/)(1) = 27m0 and 7(1) = Tm^. Substituting these val- 
ues in Eq. (|13ip . we obtain = ni — n2\f\/2 where ni 
and 712 are integers {ni = Uff, and n2 ~ 2na + for 
this case). Taking into account the arbitrariness of the 
Berry phase, vortices in G are classified by a set of two 
indices, ni and n2 mod 4. In a similar manner, we ob- 
tain /v = JT-i — 712I/I/3 for the R phase. Topologically 
distinct vortices are classified by rii and n2 mod 6. For 
C, H, I, and J, /v is described with three integers due 
to spin-gauge coupled discrete symmetries. For example, 
for the case of H, the single-valuedness condition requires 
4>{l) = 27rn0 — ATrn^/b and 7(1) = 27rn-y + 271713/5, lead- 
ing to = m — 7721/1 — "3(2 + |/|)/5. A vortex in H 
is then characterized with a set of three indices, 771, 772 
mod 2, and 773 mod 5. Note that the value of |/| in the 
above cases (phases G, R, C, H, I, and J) varies, depend- 
ing on the interaction parameters via parameters r] and 
^ (see Tables HIl and IIII| ) . Hence, the quantization unit of 
Mn/h — S{f) depends on the interatomic interaction. 



V. CONCLUSION 

We have discussed the symmetry-classification method 
based on Michel's theorem, and applied it to spin-1, 2, 
and 3 spinor Bose-Einstein condensates (BECs). We 
classify BECs having unbroken symmetries according to 
conjugacy classes of an isotropy group H, where i7 is a 
group of operations that leave the order parameter un- 
changed. For the case of spinor BECs, H \s a. subgroup of 
G = U{1) X 5*0(3) X e, where C/(l), 50(3), and 6 denote 
gauge transformations, spin rotations, and time reversal. 
For each subgroup H of G, we find an order parameter 
which is invariant under all elements of H . When H is 
large enough, the order parameter is uniquely determined 
(inert state). The obtained state is stationary regardless 
of the detailed form of the interaction energy. On the 
other hand, when H \s a small group, there exist many 
states that are invariant under H . We have characterized 
the order parameters of these states with a few parame- 
ters and found stationary points of the mean-field energy 
with respect to these parameters. The obtained order 
parameter depends on the interaction parameters in the 
mean- field energy (non- inert state). 

For spin-1 and 2 BECs, all ground-state phases are 
inert states, except for the spin-2 ncmatic (antiferromag- 
netic) phase in which two inert states and the intermedi- 
ate state between them are degenerate. For spin-3 BECs, 
there are fourteen stationary states. Among them, eleven 
states appear in the phase diagram: four of them are in- 
ert states and others are non-inert states. We have ana- 
lytically obtained the order parameters of all stationary 
states, except for the C phase, as functions of the in- 
teraction parameters. By comparing the energies of the 
obtained states, we have found a new phase (J phase) 
which exists in a very narrow region in the parameter 
space and has eluded the previous works [lol . [isj . 

Using the symmetry-classification method, we can find 
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the isotropy group of the obtained state, from which we 
see which types of topological excitations can be hosted in 
the phase. Among the obtained stationary states, A, B, 
D, and E host non-Abelian vortices since their isotropy 
groups arc non-Abelian. The mass circulation of these 
states and that of Q are quantized in units of h/{2M). 
In the other phases, the BEC has nonzero magnetization 
and the circulation of mass current is not quantized, due 
to the contribution from the Berry phase caused by spin 
textures. The difference between mass circulation and 
the Berry phase, however, is quantized as in Eq. p3ip . 

It is impossible to find a ground state that has 
no remaining symmetry, i.e., -ff = 1, by using the 
symmetry-classification method. For the case of iJ = 1, 
the symmetry-classification method amounts to solving 
Eq. ([J). For the case of spin-1, 2 and 3 BECs in the 
absence of an external field, all ground states have re- 
maining symmetries and are found by the symmetry- 
classification method. However, in the presence of an 
external field, the full symmetry G becomes smaller and 
completely broken in some phases, such as the broken- 
axisymmetry phase in a spin-1 BEC and the Zi, and 
Z3 phases in spin-3 BECs [loj . 
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Appendix A: Stationary states with discrete 
symmetries in spin-3 BECs 



The simultaneous eigenstatc of these operators is 
uniquely determined to be 

(-^,0,0,0,0,0,-^)^, (A3) 

where eigenvalues of C^z and C2X are both equal to — 1. 
The generators of the isotropy group and the mean-field 
energy for this state arc given by 

7^ = {e'^C76.,e"C2.,e^"r}, (A4) 
=c„ + gc^, (A5) 

respectively. 

2. 6*6 symmetry 

A nontrivial eigenstate of Cgz, namely, the eigenstate 
that has more than two nonzero components of the order 
parameter, is written as 

(a, 0,0, 0,0, 0,6)^, (A6) 

where a and 6 are complex numbers with |ap + |5p = 1. 
Note that we can arbitrarily choose the phases of a and h 
by applying a gauge transformation and a spin rotation 
about the z axis. Here, we choose a and h to be real 
positive numbers and write a = \J\ — r\ and h — y/rj, 
where < 77 < 1/2. Then, the mean-field energy of this 
state can be written as a function of rj as 

f (3) = 9c,(l - 2r,)2 + 12c. + 25c, _ ^^^^ ^^^^ 

The stationary point of Eq. (|A7p is rj = 1/2. Hence, 
the stationary point that has the Cg symmetry always 
possesses the Dq symmetry. 



In this appendix, we explore stationary states which 
have Dn and C„ symmetries in spin-3 BECs. 

1. Da symmetry 

The matrix representation of Cgz and C2X are given by 

Ce. -Diag[-l,e-*2V3^g-*V3^1^ 

e'''/^e*2'^/^-l], (Al) 
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0/ 



(A2) 



3. symmetry 

The matrix representation of C^z is given by 
C5z =Diag[e-'■<'^/^ e-'■'''^/^ e-'2'^/^ 

There is no simultaneous eigenstate of C^z and C2x- 

4. Cs symmetry 

In a manner similar to the case of the Cg symmetry, a 
nontrivial eigenstate of C^z can be written without loss 
of generality as 

(^2+^, 0,0,0,0,^^,0)^, (A9) 
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where —2 < 77 < 3, and the eigenvalue is 
gi4ir/5^ The order parameter in the form of 
(0, ^(3- r;)/5, 0, 0, 0, 0, ^{2 + ri)/^Y is also a nontriv- 
ial eigenstate of Csz. However, this state belongs to the 
same orbit as that of the state in Eq. (jA9p , because they 
are transformed into each other by a spin rotation and 
a gauge transformation. The mean-field energy is calcu- 
lated as a function of ?/ as 



£(3) 



^(6 + 7] -if 



which has a stationary point at 

C/3 



■n 



2(c/3 - 3c^) ■ 



(AlO) 



(All) 



The stationary point exists only when —2 < C/3/[2(c/3 — 
3c-y)] < 3. The generators of the isotropy group and the 
mean-field energy for the stationary state are given by 



^(3) _ c^(72c^ - 25c^) 



respectively. 



12(3c^ — cp) 



5. Di symmetry 



(A12) 
(A13) 



The matrix representation of C42 is given by Eq. (|118l) . 
The simultaneous eigenstate of and C2X is deter- 
mined up to an overall phase to be 



71' '''''' 7!'"^' 



(AM) 



This state has the symmetry of the octahedron, and the 
1^4 symmetry is not the largest symmetry of this state. 



6. C4 symmetry 

There are two nontrivial eigenstates of C^z- 
Case (i): The order parameter 



(^i±^,0,0,0,^^,0,0)^, (A15) 

is an eigenstate of C^z with eigenvalue i, where — 1 < 
7/ < 3. The mean-field energy of this state is written as 
a function of 77 as 

£(3)=c,r72 + g(3-ry)(7 + 77). (AI6) 



By taking derivative of Eq. (|A16|) with respect to ry, we 
find a stationary state at 



2cfi 



12c-), — cp 



(A17) 



The generators of the isotropy group and the mean-field 
energy for the stationary state are given by 



# = {e"/2C4„C2.r}, 
^(3) ^ c^(252c-y - 2hc[}) 
12(12c^ - cfs) 

respectively. 

Case (ii): The order parameter 



(A18) 
(A19) 



(0, 



1 + 11 



,0,0,0, 



2 '-'-'-v 2 



(A20) 



is an eigenstate of C^z with eigenvalue —1, where — 1 < 
77 < 1. The man- field energy of this state is given by 



2 

ll'l 1 



(A21) 



whose stationary point lies at 77 = 0. The corresponding 
state has the symmetry of the octahedron (O), and C4 is 
not the largest symmetry of this state. 



7. D3 symmetry 



The matrix representation of C^z is given by 



C^z = Diag[l,e 



i27r/3 -i27r/3 i -i27r/3 

, i, e 



(A22) 



The simultaneous eigenstate of C^z and C2X is written in 
the form of 



(, i^,0,0,e^ i±^,0,0, Ji^)-, (A23) 



where —1 < 77 < 1 and the eigenvalues of C^z and C2X 
are 1 and — 1 , respectively. The mean-field energy of this 
state is given by 

f~^'^=y(l + ^/^) + g(41 - 1877 + 417,2) 

+ ^^^^^l^(l-772)cos(2,5). (A24) 
6 



The stationary points of this function arc obtained as 

(A25) 



(b) 77 



9c^ 



5Ca + C/3 ' 



(5 = 0,7r. 



(A26) 



For case (a), the stationary state is an inert state whose 
order parameter is given by (2, 0, 0, ±iV5, 0, 0, 2)'^/3. In- 
vestigating the geometric structure of this state, we can 
find that this state has the symmetry of an octahedron. 
On the other hand, for case (b) , the the stationary state 
is a non-inert state. The order parameters at 5 = and tt 
are related to each other by a gauge transformation and 
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a spin rotation. As a result, the stationary state that has 
the D3 symmetry can be written as 



1 — 1] 



0,0, 



1 + ri 



4 ' ' ' V 2 
9c0 



0,0, 



1 — 



48cq + cp 



(A27) 
(A28) 



The generators of the isotropy group and the mean- field 
energy for this stationary state are given by 

H = {C3., e"C2., e"C/2"/V}, (A29) 



£(3) 



respectively. 



48Ca + Cj3 ' 



8. C3 symmetry 



(A30) 



There are two nontrivial cigenstates of Ca^ . 
Case (i): The order parameter 



(0, 



^,o,o,J^,o,of 



(A31) 



is the cigenstate of C^z with eigenvalue e^'^^^'^. The mean- 
field energy of this state is written in terms of i] as 



where > 0, ^ > 0, + ^ < 1, and — tt < 5 < n. Then, 
the mean-field energy can be written as a function of r/, ^ 
and 5 as 

f (3) ^co^il - 2i] + 2rf - 4e + 4^77 + 4f ) 

^(25 - 5O77 + 4177^ - 100^ + 100^?; + 100^2) 
+ 36c^e(l - e - '/) 

- 2(3ca -5cff) ^^^ - - 20 cos(2J) (A39) 
o 

which has a stationary point at (5 = and 

(lOSc^ - 12c„ - 25c/3)c^ 
9(12c2 - - c^c^) 
24CcC^ — 9CaC^ — 4c^c-y 



■n 



18(12c2 - CaC/J - CpC^) 



(A40) 
(A41) 



The generators of the isotropy group and the mean-field 
energy for the stationary state are given by 

i? = {C3.,e*"C/j/V}, (A42) 

^(3) ^ C-y(12c-yCa + l&C^Cfi - ^CgCp) ^^^^^ 
12c^ — C^Cp — CaCf) 

respectively. We have also obtained stationary points at 
5 = ±7r/2 and tt, which lie on the same orbit as that of 
the above stationary state. 



which has a stationary point at 



2(3c^ — cp) 



(A32) 



(A33) 



The generators of the isotropy group and the mean-field 
energy for the stationary state arc given by 

H = {e-*2-/3c73^^ e-*"/3[^2"^V}, (A34) 
^(3) „ c/3(72c^ - 25c/3) 



12(3c^ — c/j) ' 

respectively. 

Case (ii): The order parameter 

(a+,0,0,&,0,0,a_f , 



(A35) 



(A36) 



is the eigenstate of C32 with eigenvalue 1, where a± and b 
are complex numbers that satisfy ja+p -I- |6p -I- |a_ p = 1. 
Here, we choose a± to be real and rewrite these param- 
eters as 



b = e^'y^, 



(A37) 
(A38) 



9. D2 symmetry 

The matrix representation of C2Z is given by 

C2z = Diag[-1, 1,-1, 1, -1, 1, -1]. (A44) 

There arc two simultaneous cigenstates of C2Z and C2X ■ 
Case (i): The order parameter 



(0, 



-,0,e*' 



,0, 



1-^ n .u.ll + V , ./lzZ?,of, (A45) 



is a simultaneous eigenstate of C'2z and C2X with eigen- 
values 1 and —1, respectively. The mean- field energy is 
given by 



-f ^(l-r;2)cos(2<5). 



(A46) 



By taking derivatives with respect to rj and S, we find 
stationary states at 



(a) ry = l/4, (5 = 0, tt. 



(b) 77 



C/3 



4c/3 - 3cq 



2 



(A47) 
(A48) 



For case (a) , the stationary state is uniquely determined 
to be an inert state whose order parameter is given by 
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(0, a/3, 0, ±VTO, 0, a/3, 0)'^/4. The Majorana representa- 
tion of this state has six vertices that form a hexagon 
with the Dq symmetry. For case (b), the stationary state 
is a non-inert state. The generators of the isotropy group 
and the mean-field energy for this state are given by 



^(3) _ Cff(18Ca - 25Cff) 

3(3ca - Acfs) 
Case (ii): The order parameter 



(A49) 
(A50) 



1 - V.T 



(A51) 



is a simultaneous eigenstate of C2z and C2x with eigen- 
values — 1 for both operators. The mean- field energy of 
this state is given by 



£(3) =c„(l - 21] + 2Tf) + y|(25 - 3077 + 2677^ 



-I- sVl^cp a/77'^(1 — 77) cos 5 

+ 6(4ca - hcp)ri{l - 77) cos(2(5), 

whose stationary points are obtained as 



(a) 77 = 


1- 


0, 


(b) 7; = 






(c) 7/ = 

(d) r; = 


16 

3(4c„ - 


0, 

- 5c/3) 


8(3c„ - 




5 = 


arccos 


1-277 



(A52) 

(A53) 
(A54) 
(A55) 
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77(1 - 77) 



For cases (a)-(c), the corresponding order parameters are 
respectively given by 

(a) ^(^5, 0, 73, 0, VS, 0, Vbf, (A57) 



(b) i(V3,0,-\/5,0,-V5,0, V3)^ 



(c) ^(1,0,715,0,^15,0,1)^, 



(A58) 
(A59) 



which are all inert states and have the symmetries of Doc 
[C*") in Eq. ([TTO)) ] , O and Dq, respectively. For case (d), 
the generator of this state is given by 

H = {e'^C2.,e"C2r,e''CiyT}, (A60) 



where 



TT 4- arctan 



/I577 sin (5 



■y/1 — 7/ + \/15ri cos S 



(A61) 



Since the isotropy group generated by this H is isomor- 
phic to that generated by H in Eq. (jA49p . the order pa- 
rameter for case (d) lies on the same orbit as that of the 
state in Eq. (|A45p . 

10. C2 symmetry 



There are two nontrivial eigenstatcs of C2 ■ 
Case (i): The order parameter 



(0,a+,0,e'*V^,0,a_,0)^, 
a± = 



-^±V2' 



(A62) 
(A63) 



is the eigenstate of C2 with eigenvalue 1, where 77 > 0, 
^>0, 7;-|~^<l, and — TT < S < TT. The mean-field energy 
of this state is given by 

£(3) =l££,y(5 _ 477) + c„(l - 277 + 277^) 

-4(c„-4c^)e(l-7;-0 
2cc77(l - 77 - 2^) cos(2(5), 

which has a stationary point at 

{5cp — 6c-y){ca — 4c^) 



(A64) 



8(CaC/3 - 4C/3C^ -I- 3c^) ' 
4(CaC/3 - 4:Cj3Cj + 3c2) ' 

0. 



(A65) 



The generators of the isotropy group and the mean-field 
energy for this state arc given by 



(A56) i/ = {C2.,e"C2.r}, 



(A66) 



^(3j _ 25c„c| - 4c^c^(3c„ + 2bcp) -f 12(3c„ + Ac^y^ 



12{CaCp - ACryCp + 3c^) 



(A67) 



respectively. We have also obtained the stationary points 
at (5 = ±7r/2 and tt, which lie on the same orbit as that 
of the above solution. 

Case (ii): The order parameter in the form of 



(a,0,&,0,c,0,d)'' 



(A68) 



is the eigenstate of C2 with eigenvalue —1, where a, 6, c 
and d are complex numbers satisfying |ap + |&p + |cp + 
|c?p = 1. We have numerically minimized the mean-field 
energy of this state with respect to a, b, c and d. In the 
obtained stationary state, we can choose all components 
in the order parameter to be real and a,b,d > and 
c < 0. 
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